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Abstract 

Nuclear Magnetic Resonance (NMR) has provided a valuable experimental testbed for quantum 
information processing (QIP). Here, we briefly review the use of nuclear spins as qubits, and discuss 
the current status of NMR-QIP. Advances in the techniques available for control are described along 
with the various implementations of quantum algorithms and quantum simulations that have been 
performed using NMR. The recent application of NMR control techniques to other quantum com- 
puting systems are reviewed before concluding with a description of the efforts currently underway 
to transition to solid state NMR systems that hold promise for scalable architectures. 
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INTRODUCTION 



Nuclear spins 



"eature prominently in most condensed matter proposals for quantum com- 



puting III 0, 0, 14 , 0, y, Q, y , |9|, |10|], either directly being used as computational or storage 
qubits, or being important sources of decoherence. Fortunately, the coherent control of 
nuclear spins has a long and successful history driven in large part by the development of 
nuclear magnetic resonance (NMR) techniques in biology, chemistry, physics, and medicine 
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. The central feature of NMR that makes it amenable to quantum information pro- 
cessing (QIP) experiments is that, in general, the spin degrees of freedom are separable 
from the other degrees of freedom in the systems studied, both in the liquid and solid state. 
We can therefore describe the Hamiltonian of the spin system quite accurately; there is an 
extensive literature on methods to control nuclear spins, and the hardware to implement 
such control is quite precise. 

This readily accessible control of nuclear spins has led to liquid state NMR being used 
as a testbed for QIP, as well as to preliminary studies of potentially scalable approaches to 
QIP based on extensions of solid state NMR. The liquid state NMR QIP testbed, although 
it is not scalable, has permitted studies of control and QIP in Hilbert spaces larger than 
are presently available with other modalities, and has helped to provide concrete examples 
of QIP. Here we review what has been learnt in these initial studies and how they can be 
extended to the solid state where scalable implementations of QIP appear to be possible. 

The DiVincenzo criteria for quantum computation provide a natural starting place to 
understand why NMR is such a good testbed for QIP, and in particular, for implementing 
quantum algorithms using liquid state NMR techniques. These criteria concern (1) the 
qubits, (2) the initial state preparation, (3) the coherence times, (4) the logic gates, and (5) 
the readout mechanism. 



Quantum Bits 

Protons and neutrons are elementary particles which carry spin-1/2, meaning that in a 
magnetic field B, they have energy —fl ■ B, where the magnetic moment fl = {I is the 
spin operator) is quantized into two energy states, | |) and | t) • These states have an 
energy scale determined by the nuclear Bohr magneton fiN = eh/2mN ~ 5.1 x 10^^^ A/m^ 
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(Table HI). Since spin is inherently a discrete quantum property which exists inside a finite 
Hilbert space, spin-1/2 systems are excellent quantum bits. 





19p 


31p 


13c 


29Si 


15N 


500 MHz 


470 MHz 


202 MHz 


125 MHz 


99 MHz 


50 MHz 



TABLE I: Atoms with spin-1/2 nuclei typically used in NMR, and their energy scales, expressed 
as a resonance frequency. Frequencies are given for \B\ ~ 11.74 Tesla. 

Nuclear spins used in NMR QIP are typically the spin-1/2 nuclei of ^H, ^^C, ^^F, ^^N, 
^^P, or ^^Si atoms, but higher order spins such as spin-3/2 and 5/2 have also been experi- 
mentally investigated. In liquid-state NMR, these atoms are parts of molecules dissolved in 
a solvent, such that the system is typically extremely well approximated as being 0(10^^) 
independent molecules. Each molecule is an A^-spin system, with N magnetically distinct 
nuclei. Typically, this molecule sits in a strong, static magnetic field, Bq, oriented along the 
z axis, such that the N spins precess about z. The spins interact with each other indirectly 
via inter-atomic electrons sharing a Fermi contact interaction with two (or more) nuclei. 
The connectivity of the chemical bonds thus determines which nuclei interact. 

Since the energy scale of the interactions is weak compared with typical values of fiN\B\, 
qubits in molecules can be independently manipulated, and provide a natural tensor product 
Hilbert space structure. This structure is essential for quantum computing, and in particular, 
system scalability. 



Initial state preparation 



The energy scale of a nuclear spin in typical magnetic fields is much smaller than that 
of room temperature fiuctuations. As seen in Table U at 11.74 T the proton has an NMR 
resonance frequency of 500 MHz, whereas room temperature thermal fiuctuations are ksT ^ 
25 meV ~ 6 THz, about 10^ times larger. As the Boltzmann distribution governs the thermal 
equilibrium state of the spins 

exp 
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where Z is the partition function normahzation factor, p \ for /c^T ^ [i^\B\. Thus, 
the room temperature thermal equihbrium state is a very highly random distribution, with 
spins being in their | |) and | X) states with nearly equal probability. 

Such a highly mixed state is not ordinarily suitable for quantum computation, which 
ideally works with a system initialized to a fiducial state such as |00---0). It was the 
discovery of a set of procedures to circumvent this limitation, which made NMR quantum 
information processing feasible and interesting . The essential observation is that a 

computational procedure can be applied to p, such that the only observed signal comes from 
the net excess population in the |00 ■ • - 0) state of the thermal ensemble. One class of such 
techniques averages away the signal from all other states. This averaging can be performed 
sequentially in time using sequences of pulses which symmetrically permute undesired states, 
spatially using magnetic field gradients which prepare spins differently in different parts of 
a single sample, or by selecting a special subset of spins depending on the logical state (|0) 
or |1)) of the unselected spins. These techniques, known as temporal, spatial, and logical 
labeling, do not scale well, and only create a signal strength which decreases exponentially 
with the number of qubits realized Il6l l. 

n 

Another class of techniques is based on efficient compression [12|, and in contrast, the 
signal strength obtained is constant with increasing number of qubits realized. Indeed, only 
0(poly(n)) space and time resources are needed to initialize 0(n) spins using this method, 
which has now been experimentally demonstrated but there is a constant overhead 
factor which prevents it from being practical until n is large, or the initial temperature of 
the spins can be brought lower by several orders of magnitude. 



Coherence times 



Nuclear spins couple very weakly with the external world, primarily due to their small 
magnetic moment, and the weakness of long-range magnetic forces. Thus, typical nuclei 
in liquid-state molecules may have a Ti timescale for energy relaxation of between 1 and 
30 seconds, and a T2 timescale for phase randomization of between 0.1 and 10 seconds. 
Decoherence may occur due to the presence of quadrupolar nuclei such as ^^Cl and ^D, 
chemical shift anisotropics, fluctuating dipolar interactions, and other higher order effects. 
Though the coherence lifetimes are long, the number of gates that can be implemented is 
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limited by the relatively weak strength of the qubit-qubit couplings (typically a few hundred 
Hz at most). 

Logic gates 

In order to perform arbitrary quantum computations only a finite set of logic gates is 
required, similar to arbitrary classical computations. One such set consists of arbitrary single 
qubit rotations and the two-qubit controlled- NOT gate. We describe how each of these is 
implemented. 

The Hamiltonian describing a 2-spin system in an external field Bqz is (setting h—1) 

TL — loaIzA + ^^bIzb + 'Ha,b ■ (2) 

Here, Iz is the spin angular momentum operator in the z direction, and Ui = — 7i(l — ai)BQ, 
where 7j is the gyromagnetic ratio for spin i, which depends on the nuclear species, «j is the 
shielding constant, and 7iA,B is the spin-spin coupling. The shielding constant depends on 
the local chemical environment of the nuclei, which shields the local magnetic field, resulting 
in a shift in frequency by an amount known as the chemical shift such that spins of the 
same type (e.g. protons) can have different resonance frequencies. That spins have different 
resonance frequencies is an important requirement because it permits frequency-dependent 
addressing of single qubits. 

Spins are manipulated by applying a much smaller radio- frequency (RF) field, Bi, in the 
x-y plane to excite the spins at their resonance frequencies cui. In the rotating frame, to 
good approximation, the spin evolves under an effective field B — BiCOs{(j))x + Sisin(0)y 
(where (f) is the RF phase). The rotation angle and axis (in the transverse plane) can be 
controlled by varying 0, the magnitude of Bi and the duration of the RF. Since it is possible 
to generate arbitrary rotations about the ^-axis using combinations of rotations about the 
X and y-axis, it is possible to implement arbitrary single-qubit rotations using RF pulses. 

Two-qubit gates, such as the controUed-NOT gate require spin-spin interactions. These 
occur through two dominant mechanisms; direct dipolar coupling, and indirect through- 
bond interactions. The dipolar coupling between two spins is described by an interaction 
Hamiltonian of the form 

^^(Ia-Tb- hi a ■ n){lB ■ n)) , (3) 
5 



where h is the unit vector in the direction joining the two nuclei, and / is the magnetic 
moment vector. While dipolar interactions are rapidly averaged away in a liquid, they 
play a significant role in liquid crystal jiol . I20I 1 and solid state NMR QIP experiments. 
Through-bond electronic interactions are an indirect interaction, also known simply as the 
scalar-coupling, and take on the form 

= '^^JIa ■ Ib = 2nJhAlzB + TT J (Ia+Ib- + Ia-Ib+) , (4) 

where J is the scalar coupling constant. This interaction is often resolved in liquids. For 
heteronuclear species (such that the matrix element of the Ia+Ib- + Ia-Ib+ term is small, 
when 2t:J<^uja — ojb), the scalar coupling reduces to 

^A,_B ~ '^t^JIzaIzB- (5) 

Multiple-qubit interactions, such as the controlled-NOT (cnot) operation, may be per- 
formed by inserting waiting periods between pulses so that J-coupled evolution can oc- 
cur. For the J-coupled two-spin system, a cnot can be implemented as a controlled 
phase shift preceded and followed by rotations, given by the sequence Cab = -Rj/a(270 = 
-90)i?^Bi?^A(-90)/?^AB(180)/?j,A(90). This is shown schematically in Fig.IH 

In summary, one and two-qubit gates are implemented by applying a sequence of RF 
pulses interlaced with waiting periods. In this sense, it is perhaps interesting to note that 
the sequence of elementary instructions (pulses and delay times) are the machine language 
of the NMR quantum information processor. 



Readout mechanism 



Readout of the quantum state in NMR QIP is not the usual ideal projective von Neumann 
measurement. Instead, the system is continually read out by the weak coupling of the 
magnetic dipole moments to an external pickup coil, across the ends of which a voltage is 
produced by Faraday induction. This coil is usually the same coil as that used to produce 
the strong RF pulses which control the spins, and thus it only detects magnetization in the 
X — y plane. The induced voltage, known as the free induction decay, may be expressed as 



V{t) 



-2yntr 



e-''^'pe^^\il^^+l, 
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FIG. 1: (a) A controlled-NOT gate acting on two qubits, (b) the controlled-NOT gates imple- 
mented by a controlled phase shift gate (specified by a unitary matrix with diagonal elements 
{1, 1, 1, —1}) preceded and followed by tt/2 rotations, and (c) the pulse sequence and spin orien- 
tations corresponding to the components in (b). Note that, unlike a conventional NMR selective 
population transfer sequence, extra refocusing is required to preserve the Bell basis exchange sym- 
metry between A and B. The z rotations are implemented via x and y-rotations, which are not 
explicitly shown. 

where Ti is the Hamiltonian for the spin system, and ly operate only on the k^'^ spin, and 
Vq is a constant factor dependent on coil geometry, quality factor, and maximum magnetic 
flux from the sample volume. 

The Fourier transform of V{t) is the NMR spectrum as shown in Fig.|21for example. When 
properly calibrated, the NMR spectrum immediately reveals the logical state of qubits which 
are either |0) or Specifically, for example, if the initial state of a two-spin ^H-^^C system 
is described by a diagonal density matrix, 

"a O" 

6 

P= (7) 
c 

_0 d_ 

(where the states are 00, 01, 10, and 11, with proton on left and carbon on right and a, 
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FIG. 2: Thermal equilibrium spectrum of a 5 spin molecule. Each peak corresponds to a certain 
logical state of the remaining four spins, which is indicated by the binary numbers. The real part 
of the spectrum is shown, in arbitrary units. Frequencies are given with respect to a;j/27r, in Hertz 
21|. 

b, c, and d denote the occupation probabilities) then after a Rx{tt/2) readout pulse, the 
integrals of the two proton peaks (in the proton frequency spectrum) are given by a — c and 
b — d, and the integrals of the two carbon peaks are given by a — 6 and c — d. Both the 
proton and carbon spectra contain two peaks because of the J-coupling interaction during 
the measurement period. 

One important issue in the readout of QIP results from NMR arises because the system 
is an ensemble, rather than a single A^-spin molecule. The problem is that the output of a 
typical quantum algorithm is a random number, whose distribution gives information which 
allows the problem to be solved. However, the average value of the random variable would 
give no relevant information, and this would be the output if the quantum algorithm were 
executed without modification on an NMR quantum computer. 

This problem may be resolved {^] by appending an additional computational step to 
the quantum algorithm to eliminate or reduce the randomness of its output. For exam- 
ple, the output of Shor's algorithm is a random rational number c/r, from which classical 
post-processing is usually employed to determine a number r, which is the period of the 
modular exponentiation function under examination. However, the post-processing can just 
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as well be performed on the quantum computer itself, such that r is determined on each 
molecule separately. From r, the desired prime factors can also be found, and tested; only 
when successful does a molecule announce an answer, so that the ensemble average reveals 



the factors. Similar modifications can 
exponentially fast quantum algorithms 



)e made to enable proper functioning of all known 
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on an NMR quantum computer. 



QUANTUM CONTROL 

Implementing an algorithm on a quantum computer requires performing both unitary 
transformations and measurements. Errors in the control and the presence of noise can 
severely compromise the accuracy with which a unitary transform can be implemented. 
Quantum control techniques are used to maximize the accuracy with which such operations 
can be performed, given some model for the system's dynamics. NMR has provided valu- 
able insight into the design of schemes to control quantum systems, as the task of applying 
pulse sequences to perform operations that are selective, and also robust against experi- 
mental imperfections, has been the subject of extensive studies 25,|3l. ^ single, isolated 



quantum system will evolve unitarily under the Hamiltonian of the system. In an ensemble 
measurement (whether in space or time), an isolated system can also appear to undergo 
non-unitary dynamics, called incoherent evolution, due to a distribution of fields over the 
ensemble (27l. l28|. An open quantum system, interacting with an environment, will deco- 
here if these interactions are not controlled. The study of quantum control using NMR can 
therefore be separated into three different subsections: 

• Coherent control: how can one design RF control schemes to implement the correct 
unitary dynamics for a single, isolated quantum system ? 

• Incoherent noise: how can coherent control be extended to an ensemble, given that 
the system Hamiltonian will vary across the ensemble ? 

• Decoherent noise: how can one achieve the desired control, when coupling to the 
environment is taken into account ? 
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Coherent control 



The density matrix of a closed system evolves according to the Liouville-Von Neumann 
equation of motion: 

^ = -^[^int + Kxt,p] (8) 

where T^int is the internal Hamiltonian of the system of qubits, and Tiext the externally 
applied control fields. More specifically, extending Eqs. [2] and [4] , the internal Hamiltonian 
for a system of spin- 1/2 nuclei in a large external magnetic field Bq is : 

TV N N 

^int = Yl - afc)5o(r)/,' + 2vr ^ ^ J^jl' ■ P (9) 

k=l j>k k=l 

where —7(1 — a;fe)i?o(r) is the chemical shift of the k^^ spin. The corresponding 
experimentally-controlled RF Hamiltonian is : 

N 

Hext = Yl -lkf{r)BnF{t)e-'^^'^'^l',e^f^''>'' (10) 

k=l 

where the time- dependent functions Bjip(t) and (pit) specify the applied RF control field, 
while /(r) reflects the distribution of RF field strengths over the sample. The spatial vari- 
ation of the static and RF magnetic fields leads to incoherence j^^]. We will return to this 
in the next section. 

If the total Hamiltonian is time-independent, possibly through transformation into a suit- 
able interaction frame, the equation of motion can be integrated easily and yields a unitary 
evolution of the density matrix p{t) = f/(t)p(0)f/^(t). Given an internal Hamiltonian and 
some control resources, how can we implement a given propagator or prepare a given state? 
In QIP, it is necessary to implement the correct propagator, which requires designing gates 
that perform the desired operation regardless of the input state. We will therefore focus 
our discussion on NMR control techniques that are universal, i.e. whose performance is es- 
sentially independent of the input state, although sequences whose performances are state 
dependent can also be useful for initialization purposes. Average Hamiltonian theory is a 
powerful tool for coherent control that was initially developed for NMR [2^. Waugh and 
Harberlen have provided a theoretical framework to implement a desired effective Hamil- 
tonian evolution of a spin system over some period of time. Such a tool fits well into the 
context of QIP since it aims to implement the correct propagator over the system Hilbert 
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space while refocusing undesired qubit-qubit interactions. The basic idea is to apply a cyclic 
train of pulses P = {Pj}jLi with Il-Pj ~ ^ ^'^ system which, in its simplest form, are 
assumed to be infinitely short and equally spaced by At > 0. The net controlled evolution 
over the period T = MAt can then be expressed as 

M 

where the " toggling- frame" Hamiltonians Hk = UlHkUk are expressed in terms of the com- 
posite pulses Uk = Y['j=i Pj, k = 1, . . . , M, Uq =1. Any average Hamiltonian (up to a scalar 
multiple) can be implemented in NMR systems of distinguishable spins j30|. This work 
has been extended to correct for some experimental imperfections or uncertainties primarily 
using symmetry arguments. Composite pulses have also been used to design robust control 
sequences as they can be designed to be self-compensating for small experimental errors 



sequences as tnev can be des 
3 26, 31, ai 33, 34, 35, 3- 



Strong modulation of the spin system currently represents the state of the art in per- 
forming selective, controlled operations in large Hilbert spaces (up to 10 qubits). Strong 
modulation of the spins permits accurate selective rotations while refocusing the internal 
Hamiltonian during the RF irradiation of the spins [3^ . Figure El shows the dependence of 
the fidelities achievable on the available control resources, using this technique 0|- These 
are simulation results, assuming only unitary dynamics. The fidelities are seen to improve 
both with improved distinguishability of the spins (higher field strengths), and stronger 
modulation (increased RF power). The drawback of this technique is that it is not scalable, 
as the time necessary to find the numerical solutions grows exponentially with the number 
of qubits (see 38[ for a review of other techniques). An alternative, scalable approach that 
relies on optimization over single and pairs of nuclei has been described [39], though it does 
not appear to work as well. 

Finall y, t here is a growing body of work in the area of optimal control theory for quantum 



systems j4C 



41 



42j |. which has also used NMR as an experimental testbed. It is foresee- 



able that it might be possible to combine the ideas presented in these studies with strong 
modulation and pulse-shaping techniques to design optimal control sequences, given some 
knowledge about the system decoherence and the control parameters. 
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Maximum Fidelity Found for a Three Spin Pulse 



-•- 300 MHz 
-A 500 MHz 
■- 800 MHz 




/ A. 

y ^ " . 




/ 

■ r 
; / 







10° 10' 10^ 10^ 10"" 

Maximum Pulse Power (Hz) 



FIG. 3: Maximum fidelity achievable versus external magnetic field strengths and maximum radio- 
frequency power available for strongly modulating pulses. The frequencies shown are the proton 
Larmor frequencies at different magnetic fields (300 MHz = 7 T, 500 MHz = 11.7 T, 800 MHz = 
18.8 T). 



Incoherent noise 



Incoherent noise arises from a spatial or temporal distribution of experimental parameters 
in an ensemble measurement [27|. It is manifested in NMR in the spatial dependence 
of the Hamiltonian Ti. The density matrix at a given location in the sample still obeys 
the Liouville-Von Neumann equation where the internal and external Hamiltonian are now 
spatially dependent. Since it is the spatially averaged density matrix that is measured, the 
apparent evolution of the ensemble system is non-unitary and yields the following operator 
sum representation 43] of the superoperator 

pit)=J2P^U^it)PmKt) (12) 

i 

where Ui is a unitary operator and Pi represents the fraction of spins that experience a 
given Ui {^^Pi = 1). This incoherent evolution can be counteracted using a different set of 
techniques than those used to deal with the decoherent errors to be discussed in the next 
section. 

Hahn's pioneering work showed that inhomogeneities in the Hamiltonian could be refo- 
cused during an experiment if an external control Hamiltonian orthogonal to the inhomoge- 
neous Hamiltonian was available (see Figure In the case of an inhomogeneous but 
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FIG. 4: Decay of the transverse magnetization of a nuclear spin in a liquid state sample with large 
Bq field inhomogeneities both with and without a Carr-Purcell (CP) sequence. 



static Hamiltonian, the correlation time of the noise is infinite. This work was extended by 
Carr and Purcell to counteract long, but not infinite, correlation time noise fluctuations due 



to molecular diffusion 



45l |. Composite pulse sequences have also been used to counteract 



n 



the effects of incoherent processes j3l|, |3^ 
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M 



35 
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but often assume specific 



input states, or still need to be proven effective over the full Hilbert space of multi-qubit 
systems. Spin decoupling fits into a similar framework, and provides a means of modulating 
the system in order to average out unwanted interactions with the environment It 
las ins pire d coherent approaches to other control problems for error correction purposes 

3, HQ. 

The use of strongly modulating pulses has been extended to incorporate incoherent ef- 
fects, considering local unitary dynamics over the ensemble 0] . A priori knowledge of the 
inhomogeneity of the external Hamiltonian was used to find robust pulse sequences yielding 
a higher fidelity operation over the ensemble. This knowledge was easily obtained by spec- 
troscopic NMR techniques. Though this work focused on counteracting the main source of 
incoherent errors in an NMR experiment, i.e. RF inhomogeneity, it could easily be extended 
to compensate for experimental uncertainties such as the phase noise of the RF irradiation 
or static Bq field inhomogeneities. 
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Decoherent noise 



If the coupling between the system and the environment is weak enough and the corre- 
lation time of the noise is short, the evolution of the system is Markovian and obeys the 
following master equation: 

d\p{r)) 



dt 



tH{r) + T) \p{r)) 



(13) 



where F is the Liouville space relaxation superoperator 12]. This equation yields a non- 
unitary evolution of the density matrix so that pure states can evolve into mixed states. To 
understand and test models of decoherence, methods based on quantum process tomogra- 
phy (QPT) were developed to measure T 5^|53|, so that the dynamics of the system could 
be simulated more accurately. The full model of the system including coherent, incoher- 
ent and decoherent dynamics, has been tested extensively with a three qubit QPT of the 



Quantum Fourier Transform superoperator 



5J]. When knowledge about the noise opera- 



tors is available, quantum error correction (QEC) schemes can in principle be designed to 



allow quantum computing in the presence of imperfect contro 
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NMR has primarily been 



57l | and of fault-tolerant 



used to test the ideas of quantum error correction (QEC) 
quantum com put ations 58|. Experiments were carried out to investigate different QEC sce- 
narios , in addition to encoded operations acting on logical qubits 
Schemes to implement logical encoded quantum operations while remaining in a protected 
subspace have also been investigated [61] for a simple system made of two physical qubits 
and are still being studied for larger systems. 



QUANTUM ALGORITHMS 

Many quantum algorithms have now been implemented using liquid-state NMR tech- 
niques. The first quantum algorithms implemented with NMR were Grover's algorithm 



661 El and the Deutsch-Jozsa algorithm j68|, |69| for two qubits. The quantum counting 
algorithm was implemented soon afterwards using two-qubits A variety of implemen- 
tations of Grover's algorithm and the Deutsch-Jozsa algorithm have subsequently been per- 
formed. The two-qubit Grover search was re-implemented using a subsystem of a three qubit 



system 



7l|, demonstrating state preparation using logical labeling. A three qubit Grover 



search has been implemented, in which 28 Grover iterations were performed, involving 280 
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iwo-qubit gates 1?^. A three-qubit Deutsch-Jozsa algorithm using transition selective pulses 
73I 1 , another more advanced version using SWAP gates to avoid small couplings [tJ] , and yet 

another implementation without SWAP gates 75|. A subset of a five-qubit Deutsch-Jozsa 

algorithm has also been implemented j76|. 

The implementation of quantum algorithms reached a new level with the full imple- 



2l|. This work involved the use of 



mentation of a Shor-type algorithm using five qubits 
exponentiated permutations, combined with the quantum Fourier transform, which had been 
previously been implemented [t?!- 
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FIG. 5: The seven spin molecule used in the quantum factoring NMR experiment, showing its 
J-coupling constants, Ti and T2 relaxation times (in seconds), and chemical shifts (in Hertz) at 
11.74 Tesla. 

The most complex quantum algorithm realized to-date is the demonstration of Shor's 



algorithm using liquid-state NMR QIP methods: In this work j78|, a seven-qubit molecular 
system was used to factor the number 15 into its prime factors. This molecule, shown in 
Fig. was specially chemically synthesized to give resolvable fluorine spectra, in which the 
two ^^C nuclei, and the five ^^F nuclei, could each be addressed independently because of the 
spread of their resonant frequencies. The NMR spectra of this molecule are quite remarkable; 
for example, the thermal spectrum of ^^F spin number 1 shows 64 lines, corresponding to 
the random states of the other 6 spins (Fig. IHl). Several hundred pulses were applied, with a 
wide variety of phases, and shapes, at seven different frequencies, in this demonstration of 
the factoring algorithm (Fig. Ej). A comparison of the experimental results with numerical 
simulations suggests that decoherence was the major source of error in the experiment 
rather than errors in the unitary control, which is remarkable considering the number of 



15 



pulses applied. 




FIG. 6: Experimentally measured thermal equilibrium spectra of the NMR spectrum of fluorine 
atom number 1, in the molecule of Fig. El The real part of the spectrum is shown, in arbitrary 
units. Frequencies are given with respect to uJi/2TT^ in Hertz. 

Most recently, among NMR implementations of quantum algorithms, has been the real- 
ization of a three-qubit adiabatic quantum optimization algorithm |79|. In this work, the 
^H, ^^C, and ^^F nuclei in molecules of bromofluoromethane were used as qubits, and the so- 
lution to a combinatorial problem, maxcut, was obtained, using an optimization algorithm 
proposed by Farhi 8^ and Hogg js^. This algorithm is notable because it is fundamentally 
different in nature from Shor-type quantum algorithms, and may obtain useful speedups for 
a wide variety of optimization problems. 



Other quantum protocols 

The Greenberger-Horne-Zeilinger state and its derivative entangled states of three par- 
ticles have been studied as well. First, an effective-pure GHZ state was prepared [s^, and 
later a similar experiment was done with seven spins 3^. The claim of having created 
entangled states was later refuted based on the fact that spins at room temperature are 
too mixed to be entangled 0]. GHZ correlations have since been further studied on mixed 



states 



A quantum teleportation protocol was implemented using three qubits j85[- Superdense 
coding has been realized j86|, and an approximate quantum cloning experiment has been 
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FIG. 7: Pulse sequence for implementing Shor's algorithm to factor = 15 (for case a = 7), using 
seven qubits. The four steps through 3 correspond to different steps in Shor's algorithm. The 
tall lines represent 90° pulses selectively acting on one of the seven qubits (horizontal lines) about 
positive X (no cross), negative x (lower cross) and positive y (top cross). Note how single 90° pulses 
correspond to Hadamard gates and pairs of such pulses separated by delay times correspond to 
two-qubit gates. The smaller lines denote 180° selective pulses used for refocusing, about positive 
(darker shade) and negative x (lighter shade). Rotations about z are denoted by even smaller and 
thicker lines and were implemented with frame-rotations. Time delays are not drawn to scale. The 
vertical dashed black lines visually separate the steps of the algorithm; step (0) shows one of the 
36 temporal averaging sequences. 

implemented (an unknown quantum state cannot be perfectly copied; it can only be approx- 
imately cloned) 0]. The quantum Baker's map has also been implemented jssl . 

Several experiments have been performed in an attempt to increase the thermal polar- 
ization of nuclear spins in liquid solution, as this poses a significant challenge for scaling 
NMR quantum computers to many qubits. An algorithm approach implementin g th e basic 
building block of the Schulman-Vazirani cooling scheme has been demonstrated |l8|. High 
initial polarization of the proton and carbon in a chloroform molecule have been obtained 
by transfer from optically pumped rubidium, through hyperpolarized xenon, and a two- 
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Simulated System (S) 



Physical System (P) 



|s> $ ► |p> 

U = exp(-iHJ) Vt 

|s(T)> |p(T)> 

FIG. 8: Correspondence between the simulated and physical system. The initial state s evolves to 
s{T) under the propagator U. This process is related to the evolution of state p in the physical 
system by an invertible map (j). 

,uMt G™.ch taple„.ented on this non-the.n,a,ly polarized sy.te^ Q. In a difaant 
approach |20| para hydrogen was transformed into a suitable molecule leading to a polar- 
ization of 10% which is much larger than the thermal polarization of O(10~^). A quantum 
algorithm was subsequently performed on this molecule. Most recently, a two-spin system 
was initialized to an effective purity of 0.916 by chemically synthesizing a two-spin molecule 



using highly polarized j)ara-hydrogen 



9^. 



QUANTUM SIMULATION 

In 1982, Feynman recognized that a quantum system could efficiently be simulated by 
a computer based on the principle of quantum mechanics rather than classical mechanics 
j3| • This is perhaps one of the most important short term applications of QIP. An efficient 
quantum simulator will also enable new approaches to the study of multibody dynamics and 
provide a testbed for understanding decoherence. 

A general scheme of simulating one system by another is expressed in Figure |H1 The 
goal is to simulate the evolution of a quantum system S using a physical system P. The 
physical system is related to the simulated system via an invertible map (p , which creates 
the correspondence of states and propagators between the two systems. In particular, the 
propagator U in the system S is mapped to = (f)U(j)^^ . After the evolution of the physical 
system from state p to pt-, the inverse map brings it back to the final state s{T) of the 
simulated system. 

The first explicit experimental NMR realization of such a scheme was the simulation of a 
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FIG. 9: NMR signals demonstrate a quantum simulation of truncated harmonic oscillator. The 
solid lines are fits to theoretical expectations. Evolution of the different initial states are shown: (a) 
evolution of [0) with no oscillation (b) evolution of |0) , showing 20 oscillations (c) evolution 

of |0) + |1) + |2) + [3), showing Q. oscillation and (d) 30 oscillations. 



truncated quantum harmonic oscillator (QHO) j93|. The states of the truncated QHO were 
mapped onto a 2-qubit system as follows 
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111) 



(14) 



The propagator of the truncated QHO 



U = exp |-. Q|0)(0| + + ^|2)(2| + ^|3)(3|) ^t} (15) 

{VL is the oscillator frequency) was mapped onto the following propagator of a two-spin 
system 

Vt = exp {z (2/2 (1 + II) - 2) nr} . (16) 

Implementing this propagator on the 2-spin system simulates the truncated QHO as shown 
in Figure ini 

Quantum simulation however is not restricted to unitary dyamics. It is sometimes possible 
to engineer the noise in a system to control the decoherence behavior and simulate non- 
unitary dynamics of the system j94|. Simple models of decoherence have been shown using 
a controlled quantum environment in order to gain further understanding in decoherence 
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System 



FIG. 10: Basic model for the system, local quantum and classical environments. 



mechanisms. In one model |95|, the environment is taken to be a large number of spins 
coupled to a single system spin so that the total Hamiltonian can be expressed as 

N N 

H = uJl + ^ uj,l', + 27r ^ Ji,/i/,^ (17) 

k=2 k=2 

corresponding to the system, the environment, and the coupling between the system and 
the environment respectively (the couplings within the environment were omitted here for 
simplicity). Note that the form is identical to the weak coupling Hamiltonian of a liquid 
state NMR sample presented in the previous sections. However, the number of spins in 
a typical QIP NMR molecule is small, which makes the decoherence arising from the few 
"system-environment" couplings rather ineffective, as the recurrence time due to a small 
environment is relatively short. This can be circumvented by using a second "classical" 
environment which interacts with a small size quantum environment (see Fig. for an 
illustration of the model) 96^ . 

In this model, following the evolution of the system and the small quantum environment, 
a random phase kick was applied to the quantum environment. This has the effect of 
scrambling the system phase information stored in the environment during the coupling 
interaction and therefore emulates the loss of memory. When the kick angles are averaged 
over small angles, the decay induced by the kicks is exponential and the rate is linear in the 
number of the kicks 9^. As the kick angles are completely randomized over the interval 
from to 2-7?, a Zeno type effect is observed. Figure ^2 shows the dependence of the decay 
rate on the kick frequency: the decay rate initially increases to reach a maximum and then 
decreases, thereby illustrating the motional narrowing or decoupling limit. This 
NMR-inspired model thus provides an implementation of controlled decoherence yielding 
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FIG. 11: Simulation showing the dependence of the decay rate on the kick rate, and the onset of 
the decouphng hmit. Beyond 900 kicks/ms the decay rate decreases. 



both non-exponential and exponential decays (with some control over the decay rates), and 
can be extended to investigate other noise processes. 

A type-II quantum computer is a hybrid classical/quantum device that can potentially 
solve a class of classical computational problems [qtI . 98 |. It is essentially an array of small 
quantum information processors sharing information through classical channels. Such a lat- 
tice of parallel quantum information processors can be mapped onto a liquid state NMR 
system by mapping the lattice sites of the quantum computer onto spatial positions in the 
nuclear spin ensembles. The implementation of a type-II quantum computer using NMR 
techniques has been demonstrated in solving the diffusion equation |99|. The experimental 
results show good agreement with both the analytical solutions and numerical NMR sim- 
ulations. The spatial separation of the different lattice sites in the ensemble allows one to 
address all the lattice sites simultaneously using frequency selective RF pulses and a mag- 
netic field gradient. This yields a significant savings in time compared to schemes where the 
sites have to be addressed individually. 



CONTRIBUTIONS TO OTHER QC SYSTEMS 



NMR QIP studies have contributed significantly to enabling quantum computation with 
other physical systems. Fundamentally, this has been because of the exquisite level of 
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control achievable in NMR, which remains unrivaled. Several of these contributions are 
briefly summarized below; a complete discussion is available in the literature 38^ . 

Composite pulses: trapped Ions 

The use of composite pulses has been an important contribution of NMR to QIP. A 
single, imperfect pulse is replaced by a sequence of pulses which accomplishes the same 
operation with less error. Historically, in the art of NMR, such sequences were first invented 
to compensate for apparatus imperfections, such as frequency offsets and pulse amplitude 
miscalibrations. For example, the machine may perform rotations 

Rnie) = exp + e)n-I , (18) 

where e is an unknown, systematic pulse amplitude error. Ideally, e is zero, but in practice, 
it may vary geometrically across a sample, or slowly, with time. Using average gate fidelity 
as an error metric, this pulse can be shown to have error which grows quadratically with e. 
In comparison, consider the sequence 

BBle = R^{7r)Rs^{37f)R^{7f)R,{9) , (19) 

where R(j,{-) denotes a rotation about the axis [cos0, sin0,O]j__and the choice = 
cos^^(— 6'/47r) is made. This sequence, introduced by Wimperis |infll |. gives average gate 
fidelity error ~ 2l7r^e^/16384, which is much better than the O(e^) for the single pulse, even 
for relatively large values of e. Generalizations and extensions of this technique can help 
correct not just systematic single qubit gate errors, but also coupled gate errors 0, 13- 

NMR composite pulses have also recently been successfully employed in quantum com- 
putation with trapped ions. In an experiment with a single trapped ^°Ca ion, a sequence 
of 0(10) laser pulses was performed using a variety of phases, to implement a proper 
swap operation between the internal atomic state and the motional state of the ion, and 
a controlled-phase gate. These steps allowed the full two-qubit Deutsch-Jozsa algorithm 



to be implemented 



10l|. Composite pulses have also been used in s up er c onduct ing qubits 



demonstrating robustness against detuning in a quantronium circuit 



102| 
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Shaped pulses: superconducting qubits 



NMR also widely employs shaped pulses to achieve desired control excitations of the 
spins. Typically, this shaping is performed in the amplitude and phase domain. One goal 
of this method, for example, is to achieve narrow excitation bandwidths. To first order, the 
excited bandwidth is the Fourier transform of the temporal width of the pulse. However, 
because of the nonlinear Block response of the spins to the RF exci tatio n, the first order 
approximation rapidly breaks down for more than small tip angles 103^ . Thus, in order 
to achieve sharp excitation bandwidths for different tip angles, or uniform excitation of the 
spins over a ce rtain frequency range, a panoply of pulse shapes, such as gaussians, hermite- 



104| . and fancifully named ones, including BURP and REBURP |lOq have been 



gaussians 
designed. 

These NMR techniques are applicable to precise control of quantum systems other than 
NMR. Numerical optimization can be used to sculpt pulse shapes to provide desired unitary 



transformations 
phase qubits |l06 |. 



and shaped pulses may be useful for controlling Josephson junction 



TRANSITION TO SOLID STATE NMR 



While the liquid state studies have allowed us to explore open system dynamics and to 
develop means and metrics for obtaining control in small quantum systems, these studies 
have generally been limited to less than 10 qubits. Though the decoherence times are long 
(on the order of seconds), the strength of the spin-spin coupling (used to implement 2-qubit 
gates) is small (about 100 Hz), limiting the number of operations that can be performed. In 
addition, at room temperature the density matrix characterizing the spin system is highly 
mixed and it is necessary to use pseudo-pure states [3, ll^ • As the room temperature polar- 
ization of the sample is very small (< 1 part in 10^) , the exponential loss in signal as the size 
of the spin system grows limits the number of qubits that can be observed. Another limita- 
tion to scaling liquid state NMR techniques is the use of chemistry for frequency-dependent 
addressing. As the number of qubits increases, the number of transitions that need to be 
individually addressed grows as well. These transitions all lie within a fixed (chemistry de- 
pendent) bandwidth, making it progressively harder to address a single transition without 
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disturbing any other. While techniques such as algorithmic cooling jl?!. Il07 | and cellular 
automata 108l | schemes have been proposed to overcome some of these limitations, their 
experimental feasibility has not been demonstrated to date. 

Solid state NMR approaches allow us to obtain control over a much larger Hilbert space, 
and hold great promise for the study of many body dynamics and quantum simulations. 
The most important spin-spin interaction is the through-space dipolar coupling, which is on 
the order of tens of kilohertz in typical dielectric crystals, so that it should be possible to 
implement a large number of operations (perhaps 10^) before the spins decohere. Moreover, 
in the solid state, the spins can be highly polarized by techniques such as polarization transfer 



from electronic spins . The increased polarization allows an exploration of systems with 
a larger number of qubits, and also allows preparation of the system close to a pure state. 
While traditional solid state NMR techniques rely on chemistry for addressing, it is possible 
to introduce spatial addressing of the spins using extremely strong magnetic field gradients 
that produce distinguishable Larmor precession frequencies on the atomic scale, or via an 
auxiliary quantum system such as an electron spin, a quantum dot or even a superconducting 
qubit that is coupled to the nuclear spin system. For instance, entanglement between an 



electron and nuclear spin in an ensemble has recently been demonstrated 

A variety of architectures have been proposed for solid state NMR quantum computing, 
a few of which are enumerated below : 

1. Cory et al. proposed an ensemble solid state NMR quantum computer, using a large 



number of ra-qubit quantum processor molecules embedded in a lattice [5|. The proces- 
sors are sufficiently far apart that they only interact very weakly with each other. The 
bulk lattice is a deuterated version of the QIP molecule, with no other spins species 
present. Paramagnetic impurities in the lattice are used to dynamically polarize the 
deuterium spins, and this polarization can then be transferred to the QIP molecules 
using polarization transfer techniques. The addressing is based on the chemistry of 
the processor molecule. 

2. Yamamoto and coworkers have proposed another ensemble NMR processor, using an 



isotopic_aL 
lattice 



engineered silicon substrate containing ^^Si spin chains in an ^^Si or '^°Si 
|. The ^^Si has spin 1/2 while ^'^Si and ^''Si have spin 0. A microfabricated 



(dysprosium) ferromagnet is used to produce extremely strong magnetic field gradi- 
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ents to create a variation of the nuclear spin Larmor frequencies at the atomic scale. 
Detection is performed using magnetic resonance force-microscopy. 

3. The use of an N-V defect center in diamond coupled to a cluster of ^'^C spins as a 
quantum processor has been proposed by Wrachtrup and coworkers The hypefine 
coupling between the electron spin in the defect and the carbon nuclei allows these 
carbon nuclei to be addressed individually. Using a combination of optically detected 
electron nuclear double resonance and single molecule spectroscopy techniques, they 
suggest that it should be possible to both prepare pure spin states of the system as well 
as directly detect the result of a computation by performing a single spin measurement. 

4. Suter and Lim propose an architecture based on endohedral fuUerenes, encapsulating 
either a phosphorus or nitrogen atom, positioned on a silicon surface Q]. The Cgo 
cages form atom traps, and the decoherence time of the phosphorus or nitrogen nuclear 
spin is consequently very long. Switched magnetic field gradients that can produce ob- 
servable electron Larmor frequency shifts on nanometer length scales, in combination 
with frequency selective RF pulses are used to address the spins. Each site repre- 
sents two physical qubits, the electron spin of the fullerene and the nuclear spin of the 
trapped atom, and are used to create a single logical qubit. Two qubit gates between 
different fullerene molecules are mediated by the electron-electron dipolar coupling. 



While simple gates have been demonstrated in solid state NMR [110], the fidelity of these 
gates is low, and significant experimental challenges remain. The decoherence mechanism 
in the solid state is primarily due to the indistinguishability of the chemically equivalent 
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nuclear spins. The Hamiltonian of a homonuclear spin system is U] 



n = ujJ2il + Yl^^^ (3^^^' - ■ (20) 

i i<j 

where the first term corresponds to the Zeeman energy and second term is the truncated 
dipolar Hamiltonian in strong magnetic fields. While the total dipolar Hamiltonian com- 
mutes with the total Zeeman Hamiltonian, the Zeeman and dipolar terms do not commute on 
a term by term basis. The phase memory of the spins is scrambled as they undergo energy- 
conserving spin flips with other spins in the system. Control of the dipolar interaction is 
therefore an essential element of any NMR solid state proposal. An important example of 
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FIG. 12: (a) Free induction decay of a single crystal of calcium fluoride, at two different crystal 
orientations (b) Decay of the ^^F signal while the dipolar coupling between the nuclear spins is 
decoupled. (The sequence was modified to show an oscillation, so the decoupling is not optimal 
here.) 



the precision of control is the ability to effectively suppress all internal Hamiltonians and 
preserve quantum information. Figure a) shows the measured free induction decays for 
a crystal of calcium fluoride at two different crystal orientations, while Figure ITW b) shows 
the signal while refocusing the dipolar interaction. It is therefore possible to experimentally 
extend the coherence times of the ^^F spins in a single crystal of calcium fluoride from 100 
fis (no modulation), to 2 ms using standard NMR techniques, and finally to 500 ms using 
recently developed methods |lll| . This represents an increase by approximately 4 orders of 
magnitude. In the limit of perfect coherent control of the dipolar couplings, it should be 
possible to significantly further extend the coherence time of the spins. In addition to im- 
proving coherent control, such studies provide insight into the next important contribution 
to decoherence. 

The decay of the observed FID in Figure IT^ a) is due to the mutual dipolar couplings 
of the spins. These couplings produce correlated many spin states that are not directly 
observable using standard NMR techniques. However, using multiple quantum encoding 



techniques (ll2| |. it is possible to directly measure the growth of the spin system under 



the dipolar coupling. The truncated dipolar Hamiltonian shown above is a zero quantum 
Hamiltonian when examined in the basis of the quantizing Zeeman Hamiltonian as expected. 
However, we can requantize the system in another basis (such as the x-basis for example) 
via a similarity transformation, and explore the growth of multiple quantum coherences in 
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FIG. 13: Multiple quantum encoding, combined with evolution reversal sequences allow us to follow 
the growth of the correlated spin states during the course of a free induction decay. It is seen that 
even when the macroscopic signal appears to have decayed away, the spin system remains highly 
coherent, and states involving up to 20 correlated spins are observed. 



this new basis. The dipolar Hamiltonian in the x-basis is 

^ = E {2^^ - ^ (^^+^- + ^-^+) } - i E (^+^^+ + ^-^-) (21) 

i<j ^ ^ i<j 

and is thus seen to contain both zero and double quantum terms. It is possible to directly 
observe the growth of these x-basis coherences. Figure ^1 shows the results of this exper- 
iment, illustrating the growth in the number of the correlated spins from 1 to about 20 in 
the first 150 /is following the application of a 7r/2 pulse. 

In an early demonstration of the capability to explore many body dynamics in spin 
systems, solid state NMR techniques have been used to directly measure t he rate o f spin dif- 
fusion of Zeeman and dipolar energy in a single crystal of calcium fluoride As seen 
from the Hamiltonian above, these are both constants of the motion and are independently 
conserved. Spin diffusion is a coherent process caused by the mutual spin flips induced by 
the di pola r coupling between spins, that appears diffusive in the long-time, long- wavelength 
limit It is estimated that up to 10^*^ spins are involved at the long timescales explored 
in these experiments. It was found that while the measured diffusion coefficients for Zeeman 
order were in good agreement with theoretical predictions, the diffusion of dipolar order was 
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observed to be significantly faster than previously predicted. While these experiments were 
performed in highly mixed thermal states, future experiments planned at low temperatures, 
and high polarizations should enable a more complete exploration of the large Hilbert space 
dynamics. 

CONCLUSIONS 

NMR implementations of QIP have thus yielded a wealth of information by providing 
experimental realizations of a number of proposed schemes. This in turn has guided our 
understanding of the relevant issues involved in scaling these testbed systems up in size. The 
methodologies developed are relevant across most of the physical platforms that have been 
proposed for QIP, and manifestations of this "cross-fertilization" are beginning to appear in 
the literature. 

Sohd state NMR holds great promise for scalable QIP architectures. The efforts currently 
underway to characterize and control large spin systems are essential to determining how the 
methodologies of control and system decoherence scale as a function of the system Hilbert 
space size. 
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